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1. Introduction 

The lightest supersymmetric particle (LSP) that arises in models featuring supersymmetry 
(SUSY) at low energy with the conservation of i?-parity can be a cold dark matter (DM) 
candidate [1]. The requirement that the LSP is neutral and produced in the required 
abundance in the early Universe sets strong constraints on the minimal supersymmetric 
standard model (MSSM) [1] or its constrained version (CMSSM) [2]. In the CMSSM, the 
LSP is often the lightest neutralino with a large bino component [3, 4, 5]. In particular, 
the cases of a neutralino nearly degenerate with the lighter stau [6] or stop [7] are often of 
cosmological importance, but parameter choices for which such charged particles themselves 
are the LSP are thought to be excluded on astrophysical grounds. In fact, very strong 
bounds on the presence of electrically charged and colored relics have been obtained using 
the unsuccessful searches for exotic nuclei: strong bounds on electrically charged relics have 
been obtained up to masses of the order 10 8 GeV, while masses greater than a few TeV 
for colored thermal relics remain consistent with observation [8]. 1 

However, if the LSP does not belong to the usual MSSM spectrum, other possibilities 
arise. One such attractive scenario involves supersymmetric models implementing the 
Peccei-Quinn mechanism for solving the strong CP problem [9] . In this class of models the 

1 Note that in any case for stable charged particles not excluded by searches the overclosure bound 
applies, thus strongly limiting the possibility of such heavy thermal relics, except if some mechanism for 
suppressing its number density is at work. 
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axino, the fermionic superpartner of the axion, can naturally be the LSP and constitute the 
dominant component of DM in the form of warm [10, 11, 12, 13, 14] or cold dark matter 
(CDM) [15, 13, 16]. (See also [17, 18].) This is because, unlike for the neutralino, the mass 
of the axino is not directly determined by the soft SUSY-breaking terms and can be much 
smaller [19]. 

Axinos can be efficiently produced in the early Universe through several possible pro- 
cesses. A class of thermal production (TP) processes involves scatterings and decays of 
particles in the primordial plasma. Alternatively, in non-thermal production (NTP) the 
next-to-lightest supersymmetric particle (NLSP) decays after first freezing out from the 
plasma. In addition, one can think of other possible production mechanisms, e.g. from 
inflaton decay, but they are much more model dependent and not necessarily as efficient. 
The above processes are supposed to re-generate the axinos, after their primordial popula- 
tion has been diluted as a result of inflation and subsequent reheating with Tr < f a (with 
f a ~ 10 11 GeV being the Peccei-Quinn scale) in order to avoid overclosure. Otherwise, 
the axinos would have to be very light, m a < 0.2 keV [12], with an update in [13], thus 
constituting warm DM. 

In [15, 13, 16] we conducted a detailed investigation of the axino LSP as CDM in 
the KSVZ-type (hadronic) axion models [20] coupled with the MSSM. As summarized in 
Figs. 8 and 9 of [16], we showed that, for m a > 100 keV and Tr < 5 x 10 6 GeV, the axino can 
constitute cold DM, while at a lower mass range and larger Tr it could be a warm or even 
hot DM relic. At Tr > 10 4 GeV NTP is typically subdominant while TP processes allows 
for a narrow region of m a to satisfy the relic abundance constraint VL a h 2 ~ 1. On the other 
hand, in the region of m a > lOMeV and Tr < 10 4 GeV, the NTP mechanism of neutralino 
(or other NLSP) decay after freeze-out often plays a dominant role in producing enough 
axinos, while at larger Tr TP does this too efficiently. Contrary to the case of a gravitino, 
the other popular alternative SUSY candidate for the LSP and CDM, the lifetime of the 
NLSP is much shorter than 10 4 sec and therefore avoids the very strong bound coming from 
photo-destruction of light elements after nucleosynthesis. For this reason, only the much 
weaker constraints from hadronic showers play a role for the axino LSP scenario. However, 
if the squarks are much lighter than the gluino, their decays in TP processes can reduce 
an upper bound on Tr further [16]. 

In this paper, we turn our attention to investigating the implications on the allowed 
spectra of the CMSSM of assuming the axino as the LSP and CDM, with the right amount 
of the cosmological relic abundance and with other constraints, especially from nucleosyn- 
thesis, satisfied. This is of interest because the standard paradigm is highly constrained, 
especially at not too large values of tan (3 on which we will concentrate here. While this 
makes the CMSSM highly predictive, a question arises as to whether viable alternative 
scenarios exist which would lead to very different predictions and to possibly relaxing 
the stringent cosmological bounds on the CMSSM without abandoning a supersymmetric 
explanation for CDM. 

As we will see, very different regions of the CMSSM parameter space will often become 
allowed. In particular, it will be possible to have either an electrically neutral NLSP 
(the lightest neutralino) or a charged one (the lighter stau). The NLSP will appear in a 
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collider detector as a stable lightest ordinary super symmetric partner (LOSP). Since axino 
couplings are suppressed by the Peccei-Quinn scale, the NLSP lifetime is long compared 
to the timescales of relevance for accelerator searches. Therefore, the usual mass bounds 
arising from the decay of the sparticle into the lightest neutralino are not applicable. 
Depending of the case of a chargino or slepton L(0)SP, the bounds become a bit relaxed or 
more stringent. For the case of a "stable" charged L(0)SP, a general lower bound of about 
99.4 GeV has been obtained at LEP [21]. This may, once superpartners are discovered at 
the LHC, allow one to distinguish the case of the axino from the usually assumed case of 
the lightest neutralino as the true LSP. 

In this context it might be worth commenting on the gravitino as the LSP and CDM. 
In a recent study [22] it has been concluded that one can find limited regions of the CMSSM 
where the gravitino LSP would be allowed but that, in these regions, NLSP decays typically 
provide too little relic abundance for the gravitino to fall into the favored range, if one 
does not want to spoil nucleosynthesis predictions. However, thermal production processes 
involving gluino decay and scattering may rectify this [23]. 

In the following, we will first summarize the effective axino-neutralino couplings and 
compute effective axino-stau couplings and corresponding NLSP lifetimes in the MSSM 
coupled with the KSVZ axion models. We will then discuss bounds on the NLSP number 
density from primordial nucleosynthesis constraints and summarize the calculation of axino 
relic abundance. We will next concentrate on the CMSSM and will first illustrate the 
cosmologically allowed parameter space in the CMSSM in the standard scenario with the 
neutralino LSP. Next we will re-compute the relic abundance in the CMSSM with the 
axino as the LSP and will investigate the resulting changes in the cosmological bounds on 
the (m 1( / 2 ,rao) plane for several values of the axino mass and Tr. 

2. The Axion Multiplet Interactions 

In the Peccei-Quinn (PQ) solution to the strong CP problem, a complex scalar field is used 
to break the global U{\)pq at a high scale f a ~ 10 11 GeV. Its Goldstone boson component, 
the axion, plays the role of a dynamical field Oqcd and relaxes at the origin after the QCD 
phase transition [24, 25]. After chiral symmetry breaking, the axion acquires a tiny mass 
from instanton effects [26]. As before, in this work we concentrate on the KSVZ-type axion 
type of interactions [20]. In the DFSZ case [27] in general the mixing between the axino 
and the neutralinos can become substantial and enhance the couplings discussed here. 

In the context of supersymmetry the axion field becomes a chiral multiplet [28] which 
contains not only the pseudoscalar axion and the scalar saxion with .R-parity +1, but also 
one R = — 1 state, their fermionic superpartner, the axino a. The interaction of the axino 
with gluons and gluinos proceeds via diagrams that are the super symmetric analogues of 
the fermionic triangle anomaly diagrams for the heavy states Q. After integrating out 
the heavy KSVZ quarks Q and squarks Q, where g ~ f a , there arises an effective 
dimension-5 interaction term 



where b = 1, . . . , 8, g is the gluino and G is the strength of the gluon field, while N is the 
number of flavors of quarks with the Peccei-Quinn charge, and is equal to 1 (6) in KSVZ 
(DFSZ) models. For the remainder of this paper, we understand f a to mean f a /N. 

The coupling (2.1) and the corresponding axion coupling can also be written in a 
supersymmetric form as a Wess-Zumino term in the MSSM superpotential 

Wwz = 1^Y A Tr [W a W a ] , (2.2) 

where A is the axion chiral multiplet, W a the vector multiplet containing the gluon, and 
the trace sums over color indices. An interesting feature is that the coupling above is only 
determined by the QCD anomaly of the heavy states and is not subject to renormaliza- 
tion [29] or model dependence. 

In an analogous way, Wess-Zumino terms can arise also for other gauge interactions, 
depending on the charge of the heavy quark multiplet. At high energies, where all leptons 
can be considered massless, such interactions can be rotated into the U(l)y direction [13] 
and we are left to consider only the term 

W wz = a * < i aY j AB a B a , (2.3) 

where B a is the hypercharge vector multiplet and C a yy is a model-dependent factor [30], 
which vanishes if the heavy KSVZ quarks are electrically neutral. If they have electric 
charge eQ = —1/3, +2/3, then C a yy = 2/3, 8/3. This superpotential term leads to the ad- 
dition of the following effective dimension-5 interaction term to the low-energy Lagrangian 



^bb = ^^a^BB^ (2.4) 



Cosmological implications of the effective axino-gauge boson-gaugino operators (2.1) and 
(2.4) have been extensively studied in [13]. 

However, in addition to the above interactions, there is also an effective dimension-4 
coupling of the axino to fermions and sfermions [16] 

C a^ = S i disj ^j /R ^i P R/Ll^ (2.5) 

where tpj and tpj can De an Y of the SM fermions and their superpartners. This effective 
coupling, which arises at two-loop level in KSVZ models (and therefore at a one-loop level 
in the effective theory valid much below f a ), is the dominant channel for inducing the decay 
of charged NLSPs to axinos. 

The effective vertex in (2.5) and its effect on the axino abundance was computed in 
the case of light quarks in [16] where it was found that the dominant contribution was due 
to the logarithmically divergent part of the gluon-gluino-quark loop and was proportional 
to the gluino mass nig, 

L/R „ _^_rng , ( fa\ ( Q) 
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Figure 1: The Feynman diagrams contributing to the squark-quark-axino interaction. The thick 
dot denotes the effective gluino-gluon-axino vertex. 

Here the quark mass has been neglected. This result was obtained in the low energy 
effective theory given by the MSSM plus the interaction in (2.1), i.e. from the one-loop 
diagrams as in Fig. 1, where the Peccei-Quinn scale f a was used as a cut-off in the loop 
momentum integral. This procedure is justified by the fact that the effective theory is valid 
only below the scale and in practice amounts to absorbing our ignorance about the high 
energy theory into the parameter f a . In this way a model-independent computation of the 
effective coupling was performed. (Integrating the Renormalization Group Equation for the 
effective coupling above between f a and nig with vanishing boundary condition, gives very 
similar results.) Adopting the same strategy as in [16], here we extend the computation of 
the effective vertex to the case of the tau-stau-axino, including the effects of the tau mass. 

Note that the procedure above gives the dominant coupling of the axino to the sfermions 
as long as the mixing of the axino with the other neutralinos is negligible. In the DFSZ 
case, this is not the case and therefore the computation of the production of axinos in the 
thermal bath is more involved. For what concerns NTP instead, the picture is similar, 
apart from shorter NLSP lifetimes [31]. 

3. The Lighter Stau as the NLSP 

In order to study the charged NLSP scenario in the CMSSM, first we evaluate the effective 
couplings between the lighter stau and the axino. 

3.1 An Effective Stau Tau Axino Coupling 

As shown in Fig. 2, the stau couples to an axino and a tau via triangle diagrams analogous 
to those giving rise to the squark-quark-axino coupling. The effective axino-gauge boson- 
neutralino vertex in this diagram, represented by a thick dot, is given by (2.4). 

In the following we present the result for the loop integration. We keep only the 
dominant contribution and regulate its logarithmic divergence with the cut-off f a , in the 
spirit of [16]. 
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Figure 2: The Feynman diagrams giving the dominant contribution to the stau-tau-axino cou- 
pling. 

Consider first the loop containing a photon. This diagram gives rise to the following 
effective vertex involving an axino a with momentum ki, a tau r with momentum k 2 and 
a stau fwL with momentum k\ + k2, 

iur(k2)(Gl /R (h, k 2 )P L/R + Gl /L (k 1 ,k 2 )P R/L h 5 v~ a (k 1 ) (3.1) 

The "gauge" couplings Gl /R are as follows 




where i = 1, ... ,4 denotes the four neutralino species and Zij is the matrix that changes 
basis to the neutralino mass eigenstates. 

The "Yukawa" couplings G_r/l arise due to the higgsino content of the intermediate 
neutralino and are given by 

~ 7 ^ A 3y/2aL ( C aYY Z Vi \ ( Z l3 m T \ (m^\ f f a \ 

^R/L - 2^ g7r 2 ^ cos 2 e w )\ 2m w tan 6 W cos (3 J \ f a J & \m x J " 1 ' 
i-i l 

The analogues of these Pr/l couplings were not considered in the case of the squark- 
quark-axino coupling in [16] because the dominant diagrams there involved an internal 
gluino instead of the neutralino. 

In the calculation of G R ^ L , m T appears in the stau-tau-neutralino coupling but has 
been neglected in the propagator and in the final state mass. This is justified because the 
dominant scales in the loop and the kinematics, respectively m x and rrif, are both much 
larger than m T . Moreover, since m T /mw — 0.02 and Zu and Z^ cannot both be close to 
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1, we conclude that the effect of the tau mass is small also in the coupling and that G^ L 
can safely be neglected, as long as we assume tan/3 < 50. 

We turn now to the loop containing a Z-boson instead of a photon. The dominant 
contribution to the loop integral comes from the high energy scale f a , and is the same for 
both loops. In fact, adding a Z-boson mass term to the denominator leads to a suppression 
relative to the photon case of only a few percent, which we neglect. Therefore, the only 
difference with respect to the photon loop is in the couplings of the Z-boson to the tau, 
and to the neutralino and the axino. Note also that the G^, L couplings are suppressed in 

the same way as the corresponding and therefore will be neglected. 

In summary, the two loop diagrams together lead to the following effective t^/l — r — a 
couplings 

°i*-p& i^e-) <> - -> (x) * (4) <3 - 5 > 

x U + T^A ( ^) log ( ■ (3.6) 



tan W J \ f a 

In general, the NLSP is a linear combination of the two staus, f\ = cos Of tl + sin Of ?r. 
Therefore, the photon loop leads to the following effective f\ — r — a coupling 

m r (fci) [sm0fG L (k 1 ,k 2 )PL + cosflf G R (ki, A; 2 )Pr] 75^5(^1) (3.7) 

In the following, we consider the consequences of this effective vertex for the stau decay. 

3.2 Light Stau Decay 

Through the effective vertex described above, a fi can decay into an axino and a tau, with 
a width 

Ff^ar * (sin 2 0f\ G L \ 2 + cos 2 0f\ G R \ 2 ) . (3.8) 

In order to get a rough estimate of the decay rate without considering the whole 
neutralino sector, note that the two couplings contain the terms Zum Xi Zn and Zum Xi Zi2. 
The quantity Sj Zum Xi Zn is simply M\\ of the undiagonalized neutralino mass matrix 
in the flavor basis, which is just M\. Similarly, Y 1 iZnm Xi Zi2 is M±2, which is zero. In 
the expressions for Gl and Gr above, these sums are modified slightly by the logarithmic 
factors, but as a rough estimate we can assume that they do not modify strongly the 
cancellation and ignore the Za term in Gr. Then we have 

G L 3 - tan 2 W nnc ^ ,„ 

GR ~-- i-t^o w ^-°- 96GL (3 ' 9) 



and so 



Tf^ar ~ ^( G l) 2 ^- 0-069 cos 2 Of) 



^(G Lr . (3.10) 
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This implies that the width of the lightest stau does not depend too much on its compo- 
sition. By replacing Z\im Xi Zn with M\ in the expression for Gl, and noticing that the 
logarithm has a fairly weak dependence on the neutralino masses, we find that 



Gl ~ — 



CgYY 

2ir 2 V cos 2 Ow 
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) (l-tan 2 ^) (^j log 
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(3.11) 



This gives us the following order of magnitude estimate for the lifetime of the lighter stau 
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(3.12) 
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~(7.14sec) X Cl YY [^^) {-^-) (lOO^v)- 

So the lifetime of the light stau can be larger than 1 sec and in the case of NTP the decay 
takes place during or after nucleosynthesis. 

Note that the stau can also decay into a tau, an axino and a photon through the 
effective tree diagram in Fig. 3. Even though this decay mode is not suppressed by loop 
factors, it gives a much smaller contribution to the stau width, 



Ylat 



>OT7 



(47T) 



< (41 sec)" 1 Cl YY 




.100 GeV, 

where Fix) is the same as for the squark decay [16], 

1 13x 7x 2 



F(x) 



x 



(3.13) 
(3.14) 

and in the last estimate we have taken the maximal value for F at m x = m^. Note that 
F{x) drops fast for increasing x and so we can conclude that the three body decay can 
safely be neglected. 
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3.3 Constraints from Nucleosynthesis 

The stau lifetime is of the order of one second and therefore if a large population of NLSP 
Ti remains after freeze-out when NTP is dominant, then their decays could be in conflict 
with predictions of nucleosynthesis. In fact, the produced r immediately decays with the 
lifetime of 3 x 10~ 13 sec mostly into light mesons, and hadronic showers from this secondary 
decays can destroy light elements. Fortunately the stau lifetime is much shorter than 10 4 sec 
and so the much more stringent bounds coming from photo-destruction, which strongly 
constrained the gravitino LSP scenario [22], are automatically avoided here. 

Hadronic showers resulting from decays of non-relativistic particles of a lifetime shorter 
than 10 2 sec have been studied in [32]. Although in this case it is actually the relativistic 
tau that decays, we can still use the results of the previous analysis if we instead consider 
the stau as the initial particle. Below we will follow the procedure used before in [13]. 

Since the decay of the stau into an on-shell tau is by far the dominant channel, we 
can use the approximation 

Br(fi -^qq)~ Br(r -> qq) ~ 0.63, (3.15) 

where Br(r/fi — ► qq) is the hadronic branching ratio of the tau/stau. Using this expression, 
we can write the bound from hadronic destruction of the light elements as follows 

Br(r - qq)Y fl (T F ) < (B h Y) max (r fl ), (3.16) 

where Yf x = n^/s is the number density of f\ divided by the entropy density and the 
function (Bf l Y) max (Tf 1 ) can be read out from Fig. 12 of [32]. So, for example, for a stau 
with a mass of 100 GeV and a lifetime = 7 sec, we have the bound Yf x < 0.5 x 10~ 12 . 
Note that the bound disappears completely for lifetimes shorter than 0.04 sec and has a 
weak dependence also on the decaying particle mass [32]. 

To consider the most stringent scenario, we require that the out of equilibrium decay 
of the light staus is responsible for the production of axinos in sufficient numbers to match 
the present DM abundance. Since each stau produces one axino, we have (for non-thermal 
production) Y~* T = Y^iTp), and therefore 

Y d NTP < 1.6 (B h Y) max (r fl ) (3.17) 

If axinos produced are to make a sufficient contribution to the energy density to be cos- 
mological dark matter, they must satisfy the lower bound 

m a y a >0.34eV, (3.18) 

corresponding to £la,h 2 = 0.095 [33]. In the low Tr limit (where thermally-produced axinos 
do not make a significant contribution), when the f\ is the NLSP, axinos must therefore 
have a mass greater than the following bound in order to be the dominant part of the dark 
matter 

m a >0.2eV 1 (3.19) 

{B h Y ) max {Tf 1 ) 
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Figure 4: Lower bounds on ma for C a YY — 1 and f a — 10 11 GeV. 



If rria is lower, axinos cannot constitute the dominant component of dark matter, which 
could be then made up of other species, e.g. axions [34]. To illustrate the significance of 
the constraint (3.19), let us choose C a yy = 1 and f a = 10 11 GeV. The lower bound on 
rria for this case read off from [32] is shown in Fig. 4 for M\ = (solid line) and for the 
CMSSM case of degenerate neutralino and f\ (dashed line). 

It is important to note that the condition (3.19) on the axino mass does not guarantee 
that axinos are produced in sufficient numbers to constitute the dark matter, but merely 
that, if they are then the decays of parent staus would not disrupt nucleosynthesis. Later 
we will identify regions of the parameter space where NTP is dominant and the axino 
density is in the preferred CDM range. 

4. Relic Abundance of Axinos 

In computing the relic number density ria of axinos we will follow the procedure described 
in detail in [13, 16]. Here we merely briefly summarize its main points. 

As regards TP, the number density of axinos nj can be obtained by integrating the 
Boltzmann equation with both scatterings and decays of particles in the plasma included. 
Since n„ is well below the equilibrium one for Tr <C f a , we can neglect inverse processes. 
In terms of yield 




(4.1) 
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where s = (2n 2 /A5)g s ^T 3 is the entropy density, and normally g s * = g* in the early 
Universe. By changing variables from the cosmic time t to the temperature T, we can 
write the two solutions of the Boltzmann differential equation easily in integral form 



^0)=/^ ^^+ ■■■»"' (4.2) 

135-pw - £ dT H ' + ^jr* )> " , " j - < 43 > 

where we have considered the evolution from the reheating temperature after inflation, Tr, 
down to the present temperature To. 

In computing Y~j dec , in addition to the previously considered decays of gluinos [13] and 
(s)quarks [16], we now include for consistency the decays of the staus that are generated by 
the couplings (3.5) and (3.6), even though their effect is usually negligible. Also scatterings 
arising from these effective couplings invariably give subdominant contributions [13, 16] and 
will be neglected here. The relic abundance due to thermal production is then calculated 
by using the formula 

'nj p h 2 



m a y a ~0.36eV^-^-J . (4.4) 

In order to evaluate the axino relic abundance generated through NTP processes we 
first compute the relic number density of NLSPs after they freeze out from the plasma. 
In the case of the neutralino, we include all tree-level two-body neutralino processes of 
pair-annihilation and co-annihilation with the charginos, next-to-lightest neutralinos and 
sleptons. This is a standard method which allows us to accurate compute £l x h 2 in the 
usual case when the lightest neutralino is the LSP. We further extend the above procedure 
to the case when it is the lightest stau f\ that is the NLSP. We include all stau-stau 
annihilation and stau-neutralino co-annihilation processes. In both cases we solve the 
Boltzmann equation numerically and properly take into account resonance and new final- 
state threshold effects. The procedure has been described in detail in [35]. 

Since all the NLSPs subsequently decay into axinos, a simple relation holds 

^ntp = J^_ n (4. 5) 

"1-NLSP 

In the following, we will first illustrate the impact of the cosmological constraints on the 
plane (mi/2,mo) in the CMSSM without the axino, in which case the LSP is normally the 
lightest neutralino or the stau. 



5. DM in the CMSSM without the Axino LSP 

In contrast to the general MSSM, mass spectra of the CMSSM are tightly inter-related. 
This is because the model is defined in terms of only the usual five free parameters: tan f3, 
the common gaugino mass m^, the common scalar mass too, the common trilinear soft 
scalar coupling Aq and sgn(/i) - the sign of the supersymmetric Higgs/higgsino mass param- 
eter [i. For a fixed value of tan (3, physical masses and couplings are obtained by running 
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Figure 5: The (mi/2,mo) plane for the standard scenario without the axino LSP, and for 
tan/3 = 10, Aq = and [i > 0. In the leftmost (red) region marked "excluded by LEP" the chargino 
(to the left of a solid line) or the light Higgs is too light. In the tiny dark green band marked 
u fl x h 2 (2a)" 0.095 < Q x h 2 < 0.130, while in the orange region labeled u tt x h 2 too large" Q x h 2 
exceeds the upper limit. In the white region fl x h 2 is less than the favored 2 a range but otherwise 
allowed. 

various mass parameters, along with the gauge and Yukawa couplings, from their common 
values at Mqut down to mz by using the renormalization group (RG) equations. The large 
top quark Yukawa coupling tends to push m 2 H2 to negative values around the electroweak 
scale. Therefore electroweak symmetry is spontaneously broken and fi 2 is determined, but 
not its sign. 

Experimental and cosmological constraints on the CMSSM are usually presented in the 
( m i/2i m o) plane for various representative choices of tan/3 and other relevant parameters. 
In this paper, we will not attempt a detailed study of the CMSSM. Instead, we will aim 
at demonstrating that very different patterns arise by assuming the standard scenario and 
that with the axino LSP. 

In Fig. 5, we show the plane (^i/2; r?7 -o) Ior tan/3 = 10, Aq = and for /i > 0. In the 
absence of the axino, the lightest neutralino is the LSP above a solid red line, while in the 
wedge below it the LSP is the lighter stau f%. The red region on the left is excluded by 
LEP lower bounds on the mass of the chargino Tn^± > 104 GeV (far-most solid line) and 
the light Higgs m h > 114.4 GeV [36]. In the CMSSM without the axino LSP, the wedge 
is considered to be excluded as it would give an electrically charged stable relic. Likewise, 
most of the region of neutralino LSP is excluded due to the relic abundance ^ x /i 2 exceeding 
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tan/3 = 10, A = and fi > 0. In the left panel, T R = 200 GeV, in the right, 400 GeV. The dark 
green (orange, white) regions correspond to 0.095 < ilg/i 2 < 0.130 (f2a/i 2 too large and excluded, 
too small but otherwise allowed). 



current bounds. For a stable relic we will require that its relic abundance falls into the 2 a 
range 

0.095 < n C DMh 2 < 0.130 (5.1) 

which follows from combining WMAP results [33] with other recent measurements of the 
CMB. In Fig. 5, the constraint (5.1) is fulfilled only in a very narrow (green) band just 
above and along the red line. In this case, neutralino pair annihilation is inefficient because 
the sleptons and squarks are already too heavy to reduce r2 x /i 2 through i-channel processes. 
Nor are the heavy Higgs scalars H and A light enough to allow for resonant enhancements, 
although this becomes the case at very large tan /? > 50 which we do not consider here. The 
mechanism that comes to the rescue is the neutralino co-annihilation with f% [6], which 
however, for tan (3 = 10 considered here, is efficient only for mp2 < 500 GeV. In the narrow 
white region below the green band, ^ x /i 2 is less than the favored range of (5.1). Except 
for these two, the rest of the (m^ 2 ,mo) plane is cosmologically excluded. 

Clearly, the standard paradigm is highly constrained. We will now proceed to compare 
it with the case of axino as the LSP and CDM. As we will see, the form of the cosmological 
constraint will in general become not only very different but also will strongly depend on 
the axino mass and the reheat temperature. 

6. DM in the CMSSM with the Axino LSP 

As summarized in the Introduction, if the axino is light, < 10 MeV (while still remaining 
a cold DM candidate, i.e. rria > 100 keV), then thermal production is the dominant source 
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of axinos, and non-thermal production can be neglected. The reheat temperatures favored 
by cosmology (such as to give n- a h 2 ~ 0.1) are in the few hundred TeV range, independently 
of and too for any sensible values of these parameters. 

If the axino mass is increased to 1 GeV, non-thermal production is typically still neg- 
ligible. However, the cosmologically favored range Qah 2 then requires Tr < 0(1 TeV) at 
which point squark and slepton decays in TP start playing a major role. This can be seen 
in Figs. 7 and 9 of [16] in the case of squark decays. In particular, one can see in Fig. 9 
of [16] that, at fixed rria, as squark masses decrease, a maximum Tr allowed by 9,~ a h 2 < 0.13 
also decreases. A convenient way to present this is to invert the reasoning and to consider 
the cosmologically favored ranges of squark and slepton masses at fixed Tr. Since these 
increase with m 1 / 2 an d Too, this will be reflected in cosmologically favored regions of the 
(mi/2,mo) plane at given fixed values of to„ and Tr. 

In the left panel of Fig. 6 the (m^2,Wo) plane is shown for a reheat temperature 
of 200 GeV, while on the right this is increased to 400 GeV. In both panels mj = 1 GeV. 
Other relevant parameters have been set to the following values: f a = 10 11 GeV, tan (3 = 10, 
Aq = 0, and fx > 0. The red line divides the neutralino and stau NLSP regions. Below 
the line, the stau is the NLSP, while above it the neutralino. As in Fig. 5, the red band 
along the too axis is excluded by imposing the LEP bounds on the chargino and Higgs 
mass. The cosmologically favored range is coloured green, while the region excluded by 
the requirement that axinos should not overdose the Universe (Qah 2 < 0.130) is marked in 
orange. The white regions are cosmologically allowed, but not favored, i.e. VLah 2 < 0.095. 

In both panels, it is clear that thermally produced axinos exclude a region closer to the 
too axis, while there is a cosmologically favored strip further out. As Tr is increased, the 
excluded region grows and the favored region is pushed farther out to the right. This can 
be understood by again examining Fig. 7 of [16]: at fixed sfermion and gluino masses, as Tr 
is taken to larger values, the corresponding yield and therefore (at fixed m^) Slah 2 increase. 
The important point to note here is that a portion of the cosmologically favored region 
lies in the stau NLSP wedge, which is traditionally thought to be excluded. Note that the 
bounds of Fig. (4) are not applicable here because NTP is negligible in this region and 
also that other BBN bounds are irrelevant for cold axinos [13], similarly as for neutralino 
CDM. 

We now consider the other limit in which the axino is as heavy as possible, but slightly 
lighter than the NLSP in order to avoid a strong phase space suppression in the NLSP decay. 
In this regime non-thermal production is no longer negligible. In fact, at low enough Tr, 
as the yield due to TP becomes too small, it is NTP that starts providing enough axinos. 
Again, we refer the reader to Figs. 7-9 of [16]. We plot the case of m„ ~ tonlsp, in Fig. 7. 
In the left panel the reheat temperature is 50 GeV, while on the right this is increased to 
200 GeV. The other parameters are set to the same values as in Fig. 6. 

In the left panel, in which Tr = 50 GeV, the thermal production of axinos now plays 
a subdominant role. This is reflected in the fact that the green region of favored £lah 2 
in the neutralino NLSP section, at too ~ 100 GeV, coincides with the green region of 
the neutralino LSP in Fig. 5 as a result of the relation (4.5). For the same reason, at 
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Figure 7: The (toi/2,TOo) plane for an axino only slightly less massive than the NLSP, with 
f a = 10 11 GeV, tan/3 = 10, A = 0, and fi > 0. In the left panel, T R = 50 GeV, in the right, 
200 GeV. The notation follows that of Fig. 6. 

m i/2 — 700 GeV, the green region starts to extend further down into the stau NLSP wedge, 
along the line dividing the stau and neutralino NLSP sections, as co-annihilations between 
the neutralino and the stau keep the LSP density within the expected range. Further to 
the right, the co-annihilation ceases to be effective but stau decays provide enough f^/i 2 
for a limited range of . It is interesting to note that the majority of the cosmologically 
favored region now lies in the stau NLSP wedge. In this region the lower bound on the 
axino mass given in Fig. 4 as CMSSM applies and it is in our case satisfied. Note also, that 
even at such low temperature, the NLSP is completely thermalized in the whole (mirt, mo) 
plane. 

In the right panel, in which Tr = 200 GeV, TP can no longer be neglected, and in fact 
thermally produced axinos overdose the Universe if the neutralino is the NLSP. However, 
there is still a green region at m\/2 — 1900 GeV in the stau NLSP region, in which TP and 
NTP axinos together form the cosmologically favored abundance of CDM. 

The left panel of Fig. 6 and the right panel of Fig. 7 allow one to see the effect of 
increasing m„ at fixed Tr. The cosmologically favored band due to TP axinos moves to the 
right along the TO1/2 axis, in a similar way to the effect of increasing Tr. This amounts to 
increasing sfermion and gluino masses which in turn will increase the importance of NTP 
relative to TP, as can be seen in Fig. 7 of [16]. As a result, in the right panel of Fig. 7, NT 
production gives too much relic density in the cosmologically favored band for TP axinos, 
except for the small region in the corner of small tuq and large m-i/2- in the left part of this 
region, TP axinos provide the required abundance, and as you move to the right within the 
region TP axinos become subdominant and NTP axinos provide the necessary abundance, 
before eventually producing too many of them. 
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If one increased the reheat temperature further up one would exclude the region where 
NTP axinos can be CDM, and pushed the region where TP axinos can perform this function 
to values of larger than those plotted here. On the other hand, decreasing the reheat 
temperature below a few GeV would have the effect of suppressing the relic abundance of 
NLSPs, so that neither TP nor NTP could provide axino CDM. We will not investigate 
this possibility here. 

7. Conclusion 

The Constrained MSSM provides a popular and predictive framework for analyzing prop- 
erties of low-energy SUSY. While theoretical assumptions and experimental bounds from 
LEP and b — ► 57 provide significant constraints on the parameter space, it is the cosmo- 
logical relic density of the LSP as cold DM and the requirement of its electric neutrality 
that rule out most of the (m l ^ 2 i m o) plane and, for tan (3 < 50, allow for only a tiny region 
in the plane. 

In this paper, we have considered the cosmological bounds on the scenario in which 
the LSP is not the usual neutralino but instead the axino. Our purpose was not to provide 
an exhaustive scan of the CMSSM but to illustrate with a few examples the very different 
patterns that arise. We concentrated on two cases of a fairly light axino m„ = 1 GeV and of 
rria ~ m-NLSP where the NLSP was either the neutralino or the lighter stau. In computing 
the axino relic abundance we included both thermal and non-thermal production processes. 
Our results are a function of Tr which in all cases had to be rather low, below a few hundred 
GeV in order to give the right amount of axino CDM. We further included constraints from 
nucleosynthesis which are important for lighter axinos and staus. 

The cosmologically favored regions consistent with axino CDM relic abundance are 
in general very different from the usual scenario. Depending on m„ and/or Tr, basically 
nearly any point in the (ra^^o) plane can become cosmologically allowed. This is true 
for both the regions where the NLSP is the neutralino or the lighter stau. This significantly 
relaxes the cosmological bounds obtained in the absence of the axino. 

From the point of view of collider phenomenology, with the lifetime of 0(1 sec), the 
NLSP would appear stable in a detector. However, if the NLSP is the neutralino, then very 
likely one would find its relic abundance (calculated under the assumption that it is the true 
LSP) to significantly exceed 0.13. The NLSP could also be electrically charged (stau), or 
even colored (stop) - the possibility we have not investigated here. In conclusion, in collider 
studies of the CMSSM parameter space, both the apparently cosmologically excluded bulk 
of the (w-i/2i ?Tl o) plane where the neutralino is the NLSP and the wedge at its bottom 
where the lighter stau is the NLSP, should be given equal attention as the narrow regions 
preferred by the standard paradigm. 
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